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ABSTRACT 


Functional forms of equations were derived for predicting 10-year periodic 
height increment of forest trees from height, diameter, diameter increment, 
and habitat type. Crown ratio was considered as an additional variable for pre- 
diction, but its contribution was negligible. Coefficients of the function were 
estimated for 10 species of trees growing in 10 habitat types of northern Idaho 
and northwestern Montana. A procedure using these coefficients for calculating 
the components of current annual volume increment from diameter, height, and 
past radial increment is described, and its FORTRAN implementation listed. 


INTRODUCTION 


How fast do trees grow in height? The answer requires a great deal of detail about 
the tree's environment and about the tree itself, or some specific growth measure on the 
tree that indicates the combined effects of all environmental factors and the tree's 
individual characteristics. Diameter increment is such a measure. It is much more 
easily measured than height increment, and responds to the same growth determinants. 


Diameter increment, however, is much more responsive to effects of stand stocking 
than is height increment. Past effects of stocking on height growth relative to diam- 
eter growth are indicated by the tree's form, as measured by the height/d.b.h. ratio, 
and by its crown ratio (live crown length/total height). In this report, we seek to 
develop a prediction equation that relates height increment to concurrent diameter 
increment, tree height, diameter, crown ratio, and habitat type (Daubenmire and 
Daubenmire 1968). 


To use this equation for calculating the height increment component of current 
volume increment of trees on inventory plots, a computer subroutine is provided in 
the appendix. For more extended prognoses, the same equation for predicting height 
increment can be imbedded in a computer model of forest stand development (Stage 1973). 


For many years Forest Survey in the Intermountain Forest and Range Experiment 
Station has calculated current volume increment from changes in height computed from 
diameter increment and diameter. Wright (1961) recommended a similar procedure. Both 
approaches were derived from a single curve of total height over d.b.h. by determining 
the differences in the curve height for current d.b.h. and for d.b.h. plus diameter 
increment. Hence, these approaches assume that all trees of a particular species fol- 
low height/diameter curves that have the same slope for a given d.b.h. 


Better prediction equations should be possible where data are available from 
direct measurements of height increment. Prediction equations presented in this report 
are based on analysis of 1,165 trees felled as part of management planning inventories 
of the Kaniksu, Coeur d'Alene, St. Joe, and Lewis and Clark National Forests in northern 
Idaho and western Montana. 


DESCRIPTION OF DATA 


Ten species distributed among 10 habitat types are represented in the data 
(table 1). Distribution of samples with respect to height and d.b.h. are shown by 
species in figure 1, and by habitat types in figure 2. Periodic diameter growth 
(inside bark), AD, was measured concurrently with periodic height growth, AH, for the 
most recent 10-year pene’ Values for height (H) and diameter at breast hei ght (p)1/ 


were determined as t 
sion estimate is required to predict the change in height as a function of height, 


diameter, and projected diameter increment (inside bark). Of necessity, the crown ratio 


used was observed at the end of the growth period. Fortunately, crown ratio changes 
very slowly within 10 years. 


The trees were measured during the inventories of 1969-1972. Consequently, their 
most recent 10-year growth is influenced by the climate and other temporal events from 


about 1959 to 1972. 


1/ Bark increment was not deducted when past d.b.h. was calculated. The slight 
overestimate of past d.b.h. can cause only a slight bias in the estiamte of height 


growth. 


ey existed at the start of the increment period because the regres- 


Table 1.--Distribution of sample trees by spectes and habitat type 
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DERIVING THE FUNCTIONAL FORM 
FOR THE PREDICTION EQUATIONS 


A set of height/diameter curves for three even-aged stands of grand fir is shown 
in figure 3. These curves are plotted solutions of a nonlinear equation for predicting 
tree height from diameter, site quality, and age of the stand. Analysis of these data 
confirmed that the allometric relation: 


In(#) = a + b 1n(D) (1) 


does indeed apply to the development of trees as age increases, but that the coefficients 
a and b depend upon the competitive status of the trees. Competitive status can be 
defined by the distance of the (#,D) point for a particular tree from the curve of 
relation (1) for dominant trees (Perkal and Battek 1955). The dashed line shows the 
trend followed by the average for the dominant and codominant trees in the stand. 

Trees in lower crown classes move more steeply upward. Conversely, dominants move 

along lines of lower slope. 
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Figure 3. Sequence of hetght/dtameter curves for even-aged stands of grand fir. 


Instantaneous growth rates can be obtained by taking the differential of 
equation (1): 


H 
oH = b D 3D (2) 
or, in its logarithmic form: 
In(dH#) = 1n(9D) + 1n(#) - In(D) + In(d) (3) 


On the other hand, periodic growth rates expressed as a finite difference derived 
from (1) would be: 


In(#j) - 1n(#,) = b(In(D2) - In(Dj)) 


which is the same as: 


In(Ho/H,) = b 1n(Do/D,) (4) 


where the subscripts 1 and 2 designate measured values at the start and end of growth 
period, respectively. 


Models were derived from (2), (3), and (4) by using various transformations of 
height, the ratio of height to diameter, and crown ratio to estimate the b parameter 
for each species. The best transformations of these independent variables were selected 
by combining them in groups of sets for screening overall combinations of the groups 
by using one set from each group (Grosenbaugh 1967). Coefficients in these alternative 
models were estimated by least-squares regression. Goodness-of-fit indices (Furnival 
1961) for the several transformations of the dependent variable were compared using the 
best regression for each transformation of the dependent variable. 


For the screening of these alternative models, data from 909 trees of 10 species 
in the northern Idaho forests were used. 


Differential Model 


The dependent variable in this model is the 10-year periodic height increment, AZ. 
Some coefficients (ej) in the model may be different for each species; other coeffi- 
cients may be constant for all 10 species. Those that vary with species have an addi- 
tional subscript 7. 


H AD AD D 
AH = — ADle,, + @2, In(CR) + c3 ate ey (5)? +25. Ini(@) + C6, FAD! 


D 


The expression in the brackets represents the estimate of b in equation (2). In compu- 
tation, C6, was the set of constant terms in the regression model. 


Table 2 is an analysis of variance showing the improvement in the regression sum- 
of-squares for successively more complex collections of variables. In this table, for 
example, comparison level 2 indicates that when either the coefficients of HAD/D or 
the constant terms are allowed to vary by species, then the fit is improved by a sig- 
nificant amount. Of these two alternatives, the former has a slight advantage. At 
comparison level 3, varying constant terms by species is of little value so long as the 
coefficients of HAD/D depend on species. However, the variable In(H) is shown to be 
needed to represent the effect of changing height on the b coefficient. At comparison 
level 4, crown ratio is still of little use, but addition of AD/D and its square 
results in a considerable improvement. 


Table 2.--Analysts of improvement of regresston model attributable to adding vartables. 
Dependent vartable = AH 


Proportion of ° 
>: remaining variance : Marginal > Mean 


Comparison : 
level - Source: |: distic: explained!/ > sum of squares]/ : square ‘ F 
1 1 1 49469 .808 
HAD/D 1 2263.450 2263.450 263.0 
2 [ #AD/D 5 5 | 9 0.065 776.592 86.288 10.0 
Lisp 9 060 728.725 80.969 9.4 
3 lisp 9 .0168 260.549 28.950 3.4 
[ in (Z) ] 1 e2Si: 1491.782 1491.782 17323 
4 In (CR) 10 048 479.677 47.967 5.6 
lisp 9 .0158 219.627 23.403 Zw 
[ ap/D, (ap/D)? | 2 .108 905..914 452.957 5Z).0 
Error 878 7556.068 8.606 
2) 


— Increase in explanatory power for regression due to adding variables to the 
model composed of those variables bracketed in the levels above. 


Finite Difference Model 


The dependent variable in this model is 1n(H,/H,). When an expression for b 
analogous to its form in the differential model is inserted in the finite difference 
model (4), we obtain: 


In(H#5/H,) = 1n(D>/Dj,) [e1, Te CP 1n(CR) tai C9 In(D5/D,) a(n (1n(D5/D,))? 


+ ¢5, In(H) + C6 /\n(D2/D,) | 


To compare the utility of the two alternative forms of the dependent variable, we 
use the maximum likelihood index of fit (Furnival 1961). The basis for the comparison 
is the standard error of estimate of the untransformed dependent variable. Standard 
errors of estimate for other transformations of the dependent variable are converted 
for comparison by multiplying them by the inverse of the geometric mean of the deriva- 
tive of the transformation. The derivative of ln (Hj/H,) with respect to AZ is: 


9 In(H>/Hy) 9 In(1. + AH/H}) 


f' (AH) = = 
JAH QAH 
Fees el 1 
Ho Ay Ho 


Accordingly, the inverse of the geometric mean is given by: 


Sat art 
[f' (4#)] 


The index-of-fit is 3.54 for the best collection of variables (through level 4) in the 
finite difference model. This index is compared to the standard error of estimate of 
the differential model, which is +2.93 feet. Hence, the differential model is superior. 


= exp[Z n/a 


Table 3 is an analysis of variance for the finite difference model with comparisons 
Similar to those in table 2. The variables in this model that are analogous to those in 
the differential model are effective as predictors to about the same degree. 


Table 3.--Analysts of improvement of regression model attributable to adding variables. 
Dependent variable = 1n(A>/#}) 


: Proportion of : 
Comparison : : > remaining variance : Marginal > Mean 


level ~ Source > d.f..: explained!/ > sum of squares!/ = Square : F 
1 1 Z 
In (D2/D,) 1 3.1589 3.1589 1707.5 
2 [ in(2/P1) ;.,., | 9 0.091 s5572 .0375 20.3 
ae 9 .031 .1365 .0152 SoZ 
tsp 
3 i. 9 .015 .0754 .00838 4.5 
isp 
[ 1n(Z) ] q BPS 7702 =f 102 416.3 
4 [ ap/p, (ap/0)2 | 2 .258 5921 2961 160.0 
5 ioe 9 .018 .04633 200515 2.8 
tsp 
Error 886 .00185 


erie in explanatory power for regression due to adding variables to the 
model composed of those variables bracketed in the levels above. 


Differential Model-- Logarithmic Form 


When the differences between the observed and predicted values of AH are plotted 
over the predicted values, the variance of the difference increases with larger values 
of the prediction. As a consequence, faster growing trees are given undue weight in 
estimating values of the coefficients. The logarithmic form of the model (3) should 
decrease the variance of the residuals associated with large predictions. The pre- 
ceding analyses showed that the allometric coefficient (b) varied with species and the 


11 


relation of height, diameter, and diameter increment. These effects could be incorpo- 
rated quite readily by multiplying each independent variable by a coefficient to be 
estimated, and by introducing a constant term for each species. In the logarithmic 
form, the possibility of zero values for AD must be considered. Measurements of diam- 
eter growth are usually recorded only to the nearest 1/20 inch. By shifting the entire 
scale up 1/20 inch, zero values can be accommodated without distorting the overall 
relation. 


Thus, the logarithmic form of the differential model is: 


In(AH) = ec, _ + ep In(AD + .05) + G3 In@) + cy, In) (5) 
L 


For the 909 trees from the northern Idaho forests, the following coefficients for 
this model gave a better index of fit (+2.73 feet) than the differential model with the 
loss of fewer degrees of freedom: 


3.4448 Western white pine 
3.1258 Western larch 
Saal Douglas-fir 
529759 Grand fir 
@1,= 3.1508 Western hemlock 
w 2.9357 Western redcedar 
3.0706 Lodgepole pine 
ye AOA Engelmann spruce 
2. 83595 Subalpine fir 
3.4380 Ponderosa pine 
Cn => 40537401 In(AD + 0.05) 
C3 = -0.29805 In(D) 
Cr =~ --OL 13170 1n(Z) 
ee O75 


When the same model was applied to the 265 lodgepole pine and Douglas-fir trees 
from the Lewis and Clark National Forest, the values of the coefficients were quite 
different from those above. The estimated residual variance for these data was 
s* = 0.1630. The difference between these two populations is apparent in the following 
averages: 


Northern Idaho Forests Lewts and Clark N.F. 


AH 6.5 ; 1.4 
AD 8 0.44 

Differences in growth rates correspond to marked differences in habitat for those 
species that are common to the two geographic sources of data. As a consequence of the 
differing habitats, the Lewis and Clark trees would show a much lower average site index 
than the northern Idaho trees. Furthermore, the Lewis and Clark trees are generally 
older. To accommodate these differences, data from the two sources were merged and the 
differential model modified to permit coefficients to vary with habitat as well as 
species. If all four coefficients were unique with respect to species and habitats, 
then there would be 204 coefficients to be estimated by a separate regression solution 
of equation (5) for each of the 51 entries in table 1. Instead, the four coefficients 
in the logarithmic form of the differential model were varied according to species and 
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habitat in repeated regression problems (table 4). The final line represents the best 
model of the series. Coefficients derived from this solution are given in table 5. 
The ratio of the mean-square deviation from regression to the mean-square deviation 
about the mean is 0.2722 for this model. That is, the model accounts for about 73 
percent of the variance of logarithm of height increment among the 1,165 trees in the 
sample. Converting the residual error to the scale of feet of height increment 
(Furnival's index), the mean-square error of estimate would be +1.97 feet. 


To illustrate the implications of this functional form for height increment in 
the context of stand development, the successive height/diameter curves were plotted 
for a typical even-aged stand of grand fir carried through four decades of simulated 
growth (fig. 4). The curves in figure 4 show that the prediction functions for height 
increment derived in this report can generate successive curves that conform to the 
general shape and level of typical height/diameter curves in even-aged stands shown 
in figure 3. Conformation of these curves for simulated stand development also depends 
on the way diameter increment and mortality change among diameter classes within the 
stand. 


Table 4.--Effect of varying coefficients by spectes and habitat on the mean square 
restdual of the logarithmic form of the differenttal model 


Number of : Variable : Mean 
coefficients : : : : : : square 
in model : uf : 1n(AD+.05) a in 2D) : In (4) : :residual 

22 H < 6 H 0.2045 

Su S e H H .2029 

32 S+H x - H .2008 

Ky S+H H ze ii .1958 

40 S H H H SIS 

41 S+H H H e 211853 


Key: * indicates a single coefficient for all species and habitats. 
H indicates a unique coefficient of the indicated variable for each habitat. 
S indicates a unique coefficient of the indicated variable for each species. 
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Table 5.--Coeffictents for estimating the logarithm of 10-year pertodic hetght growth 


Habitat 
S 
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hemlock SA0256) Sr Zolsin Sacoo 
Engelmann 
spruce 2.9080 Se 225i) 2.3509 1.6339 
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Vartable coeffictents 
Variable 
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Figure 4.--Projected hetght/diameter curves for grand fir on Abies grandis/Pachtstima 
habitat. 


Removing Bias in Logarithmic Form 


When the variable to be predicted is analyzed by making a logarithmic transforma- 
tion, a bias is introduced when the inverse transformation is used to convert the 
estimate to the original natural scale. Effect of the bias would be cumulative when 
these prediction equations are used repetitively to simulate growth of trees and stands 
through successive periods of time. By removing the bias, repetitive application will 
produce the same height from the accumulated predictions of increment as would be the 
result of cumulating observed values of the log-normally distributed height increment. 


Size of the bias depends on the residual variance, which is listed as 0.1853 on 
the last line of table 4, for the logarithmic model. 


Bradu and Mundlak (1970) devised a correction for the bias that recognizes the 
effect of sampling error on the parameters of the regression. Unfortunately, their 
method requires the values of the inverse covariance matrix which is a 41 X 41 array. 
As an alternative, a simple correction that is conditional on the sample values of the 
parameters seems adequate. To determine the corrected estimate of the mean on the 
original scale, the value of one-half the residual variance is added to the estimate on 
the logarithmic scale before making the transformation back to the natural scale. In 
this case, the amount to be added is 0.1853 + 2 = 0.09265. The effect is to multiply 
the uncorrected estimate of height increment by 1.0971. 
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APPLICATIONS 


Calculating Current Annual Increment 


Current annual volume increment is commonly calculated by measuring or estimating 
the two major dimensions of volume increment: change in tree d.b.h. (inches) and change 
in tree height (feet). Then, these changes are added to present d.b-h. and total 
height, and the volume change is determined from the difference between the volumes 
obtained from a volume table or volume equation that relates tree volume to d.b.h. and 
total height. The prediction equations developed in this report can be used to calculate 
the 10-year periodic height increment from the corresponding 10-year periodic diameter 
increment, initial d.b.h., and initial height. When past periodic diameter increment 
is used to determine future diameter increment, a problem arises because diameter growth | 
generally declines with age. Consequently, the current annual rate would differ from 
the periodic rate divided by period length. To remedy this problem, we can invoke the 
observation that whereas diameter increment declines with age, basal area increment 
usually remains constant (Spurr 1952, p. 214; Smith 1962, De, 59): 


The sequence of calculations is: 


1. Determine past increase in diameter squared. If increment was measured out- 
side bark, convert to increment inside bark. If measured for a period of y years, 
convert to a 10-year basis. 


= 72 2 * 
DDS. = (Gs = 45.) * WO/y 


2. Use past DDS to estimate future AD (10-year periodic increment): 


AD = vd.” + DDS - d, (inches) 
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3. Solve equation (5) using the set of coefficients from the column of table 5 
for the applicable habitat type, and using the first coefficient for the appropriate 
species: 


. ip{Ad) =e; + .c>) InfAD’ + 0:05) + ¢3 In(0 ;) => (Eth In (@ ;) 


AH = exp(1n(AH) + 0.09265) (feet) 


4. Convert growth estimates to diameter and height 1 year hence: 


D.. = YD,” + BKR* DDS/10 


ny 
" 


H. + AH/10 
Lt 


where BKR = (D,/d;)* is the square of the ratio of d.b.h. to diameter inside bark at 
breast height. 


A computer subroutine prepared in FORTRAN IV to carry out this sequence of calcula- 
tions is given in the appendix. 


If growth for an arbitrary interval (y) greater than 1 year is desired, then the 
division by 10 in the last two equations would be replaced by multiplication by the 
value of y/10. That is, 


- = - 

Digy VD ,2 + BKR*DDS*y/10 
= * 

Heyy = Hy + OH*y/10 


Le, 
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APPENDIX 


FORTRAN IV Subroutine 


QAO @iG@y Ge) 


DEG e GGG) Gn Gal! Gd Ga) GaGa), 


GVOlG: GQ Gi @iGiGy GiGi ty G 


DATA ASPE /-0.260831,-0- 531522 +-0.536742,-0.423059, 


1 


i 


1 


l 
1 


1 


ISP! [S SPECTES “SUBSCRIPT 
IHAB TS HABITAT SUBSCRi PT | 
D IS DIAMETER 

HLS HEIGHT 

DD IS LO-YFAR PAST DbH GROWTH 
D2 TS THE NEW DIAMETER 


H2 IS THE NEw HEIGHT 
NIMENSTON ASPE (LU) 2AdAB(iO) »BHAB{10) »>CHAB(10), BKRI 10) | 


| 

SUBROUTINE HTGT(1SPeITHABsD0,H,DD,02-H2) 
SPECIES SUBSCRIPT DEFINITION | 
wilT& PINE | 
WE STERN LARCH 

DUUGLAS FIR 1 
vRAND FIR ir 
WESTEKN HEMLOCK } 
CEDAR | 
LODGEPOLE ! 
SPRUCE | 
SUBALPINE FIR 
10 PONUDEKOSA PINE | 


ODNAON FS WN & 


—0.5692649-06 7622 12 1-00 485439 4-0.625377 ,-02717423, -02185745/ | 


HABITAT SUBSCRIPT DEFINITION 
PSEUDOTSUGA / PHYSOCARPUS | 
PSEUUVOTSUGA / SYMPHORICARPOS | 
PSEULOTSUGA / CALMAGROST IS I. 
ABIES GRANDIS / PACHISTIMA 
ABIES LASIOCARPA / PACHIST IMA 
ABLES LASICCARPA / XEROPHYLLUM 
ABLES LASIOCARPA / MENZIESIA 
ABLES LASICCARPA / VACCINIUM | 
THUJA / PACHISTIMA 
10 TSUGA HETEROPHYLLA / PACHISTIMA | 

DATA AHAB/3228270 022962 02 92053 339 94201652 44.4574], 

22976259 2625931 930 34406 13259483 93285106/ 
DATA BHAB/ .11127472 42 9658621 5261883491 90449199914 «61660659, 

© 0559213 53 924821 9198 96 15457698 922709992 59 04090 7645/ 
DATA CHAB/.002 7789974 »—O. 31462085 9-0 565384035 —0 035620483» 

— 06472275145 6 049378771 9021472478 06 37078756, —0 223159795 

-0.22855937/ 
DATA DSP/-0.19554013/ I. 
DATA BKR /10* 1.09/ | 
D2 = SQRT(ID#*2 + DD * (2. * D —- DDI * BKRIISP) / 10.) 
H2 = H + EXP(ASPE(ISP) + AHAB(IHAB) + BHAB(IHAB) * ALOG(DD + 0.05) 

+ CHAB{IHAB) * ALUS{D) + DSP * ALOG(H)) * 0.109698 
RE TURN | 
END 


OON AUF WN 
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